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Abstract 

In our recent paper, we have obtained a model solution to the problem of radiation-inefficient accretion flows 
(RIAFs) in a global magnetic field (so called, resistive RIAF model), which is asymptotically exact in outer regions 
of such flows forming accretion disks. When extrapolated inwardly, the model predicts a local enhancement of 
the vertical Poynting flux within a small radius that may be regarded as the disk inner-edge. This fact has been 
interpreted as the origin of power source for the astrophysical jets observationally well-known to be ejected from 
this type of accretion disks. Since the accuracy of the solution may become rather poor in such inner regions, 
however, the ground of this assertion may not seem to be so firm. In the present paper, we develop a sophisticated 
discussion for the appearance of jet-driving circumstances, based on a much more firm ground by deriving a global 
solution in the same situation. Although the new solution still has an approximate nature, it becomes exact in the 
limits not only of large radius but also of small radius. The analytic results clarify that the electrodynamic power 
is gathered by the Poynting flux, from outer main-disk region to feed the innermost part of an accretion disk. The 
injected power largely exceeds the local supply of work by the fluid motion. 
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1. Introduction 

The series of states of accretion disks called the radia- 
tively inefficient accretion flows (RIAF) forms an optically- 
thin, under-luminous (usually radiating at a small fraction of 
the Eddigton luminosity of each central object) branch in the 
accretion-rate vs. surface-density diagram. Another sepa- 
rate branch exists in a more optically-thick (i.e., large surface- 
density) domain and continues from the standard-disk series to 
the slim-disk series, via a thermally unstable part, as the accre- 
tion rate increases (e.g., Kato et al. 2008). Specifically for the 
RIAF theories, a more detailed description can be found, e.g., 
in Narayan & McClintock (2008). 

The main efforts to take the effects of ordered magnetic 
fields into account in the accretion disk theories may be di- 
vided into two classes. In one class, the presence in the disk 
of only a toroidal field with even polarity (i.e., the polarity is 
the same on both upper and lower sides of the equatorial plane) 
is taken seriously. The resulting magnetic pressure is added 
to gas pressure to support the disk against the vertical com- 
ponent of gravity. Further, if the a-prescription (Kato et al. 
2008) with respect to the total pressure is adopted for a vis- 
cosity, an additional viscous extraction of angular momentum 
passing through the disk plane becomes possible. For example, 
the modifications of the standard-disk theory (e.g., Begelman 
& Pringle 2007) and of RIAFs (e.g., Oda et al. 2012) have 
been discussed, respectively, in relation to some controversial 
spectral features seen in cataclysmic variables and to the state 
transitions seen in Galactic black-hole X-ray binaries. 

In the other class, on the contrary, the presence of both 
poloidal and toroidal components of an ordered field are taken 
seriously. Since the toroidal component is considered to appear 



as a result of dragging of the vertical field lines by the rotational 
motion of the disk, its polarity reverses on both sides of the 
equatorial plane (i.e., an odd polarity). Thus, the toroidal com- 
ponent develops mainly outside the disk and vertically com- 
presses the disk against gas pressure. Moreover, such a con- 
figuration makes it possible to vertically extract the angular 
momentum by the Maxwell stress. This point is essential in 
relation to the production of astrophysical jets (e.g., Blandford 
& Payne 1982; Pudritz & Norman 1983; Lubow et al. 1994) 
often observed to emanate from the vicinity of the disk inner 
edge. 

In most of the analytic models that is addressed to the for- 
mation of jets, however, the magnetic field is not treated self- 
consistently with the fluid disk. Self-consistent inclusion of an 
ordered magnetic field into RIAF states has been performed 
in a series of works by the present author (for a review, see 
Kaburaki 2012; hereafter referred to as Paper I). In this model, 
a twisted magnetic field works to extract angular momentum 
from the disk plasma, and the resistive dissipation converts the 
available gravitational energy into heat. This makes a good 
contrast with the usual RIAF models, in which only turbu- 
lent magnetic fields are included, and the fluid viscosity plays 
an essential role in converting energy and extracting angular 
momentum. Therefore, we call the former the resistive-RIAF 
model, distinguished from the latter, the viscous-RIAF model. 

It should be mentioned also that there is another series of 
studies in which the presence of an ordered magnetic field 
is treated self-consistently (for a review, see Ferreira 2008). 
Although its relation to RIAFs is not so clear, Ferreira and 
his coworkers discuss an inner region of the accretion disk 
threaded by a global magnetic field. Their main interest is in 
the local (i.e., at a given radius) mechanisms to launch magne- 
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tohydrodynamic (MHD) jets, and the details of vertical trans- 
port of energy and angular momentum are investigated. On the 
other hand, the present concern of the resistive-RAIF model is 
to show how the energy can be supplied to the jet launching 
site from wide aria of an accretion disk. 

This paper is a direct follow-up of Paper I that has been 
devoted to discuss the appearance of the Poynting flux near 
the inner edge of a resistive-RIAF, which may lead to the jet 
launching. However, the discussion was based on the inward 
extrapolation of an outer asymptotic solution whose accuracy 
is not necessarily guaranteed in the inner region. Moreover, 
the outer solution has been derived by assuming a specific con- 
dition, which we call hereafter the extended RIAF condition 
(equation [9] in Paper I or [22] below). This condition may 
seem rather arbitrary or artificial. Therefore, we give it up in 
the present paper. Instead, according to the spirit of this con- 
dition, we first obtain several asymptotic solutions in the outer 
region of an accretion disk, which are equivalent to each other 
within the accuracy to the first order in the smallness parameter 
€q (the definition will be given in the next section). 

Under the above situation, the criterion to sift a specific so- 
lution from others would be the wideness of its applicability 
range. Thus, we are naturally led to examine the behavior 
of these outer solutions in the opposite limit of small radius, 
and find that only one case among them becomes exact also 
in this limit. Namely, the selected one becomes accurate not 
only in the limit of large radius but also in that of small radius. 
Therefore, it may be called a global solution, although it is still 
an approximate one at middle radii. This finding is indeed a 
great improvement since we can discuss global operation of 
such accretion flows based on this much secure ground than 
before. Another advantage of this improved solution is that 
the expressions for all relevant physical quantities are written 
analytically in closed forms. 

The organization of this paper is as follows. In section 2, the 
variable-separated version of the governing equations are sum- 
marized, extracted from Paper I for convenience. As plausible 
examples of asymptotic solutions at large radii, four possibili- 
ties are derived in section 3 without employing the extended 
RIAF condition. By examining the behavior of these outer 
asymptotic solutions in the opposite limit of small radius, we 
find in section 4 that there is one and only one case in which 
the same expressions become asymptotically exact also in this 
limit. Full expressions for the relevant quantities within this 
global solution are derived also in this section. Using these ex- 
pressions, we calculate and discuss in section 5 the local energy 
budgets of a few types. As summarized in the final section, the 
obtained results clearly show the appearance of preferable cir- 
cumstances for jet launching. 

2. Fundamental Equations 

As for notation, we completely follow Paper I, and hence 
adopt spherical polar coordinates, (r, 8, tp). The normal- 
ized versions of the radius and co-latitude are, £ = r/r- m , and 
rj = (8 — 7r/2)/A, respectively. Here, r- m denotes the inner- 
edge radius of the accretion disk. All physical quantities have 
been expressed in the variable-separated forms, equations (16)- 



(30) in Paper 1 1 , within the geometrically-thin disk approxima- 
tion where the half opening-angle of the accretion disk is very 
small, i.e., A <^ 1. The fundamental equations for the radial- 
part functions in quasi-stationary problems (for the definition, 
see Paper I) are summarized below. 

leading equations 

• magnetic flux conservation: 



• Ampere's law: 
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• Ohm's law: 
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mass continuity: 
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equation of motion 
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D = 



(11) 

(12) 
(13) 

(14) 



1 There is a typographic error in these equations. On the right-hand side of 
equation (21), tanhr; should read tanh 2 rj. 
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• equation of state: 

P = KpT, 

subsidiary equations 

• Faraday's law: 

1 db r _ ^E v 
c dt r 



(15) 



1 dbg _ Id 
c dt r dr 

1 db v _ 1 f d 
c dt r \ dr 

• charge density: 

~ = Ee_ 
Anr 



(rEg) 



■2E, 



(16) 
(17) 
(18) 

(19) 



One of the characteristic aspects of our treatment is that 
Ohm's law is used directly, without substituted into Faraday's 
law. The poloidal and toroidal magnetic-Reynolds numbers, 
equations (8) and (9), have been introduced in rewriting Ohm's 
law. The final expression for 5R p can be derived by using equa- 
tion (73) of Paper I. 

The expression (11) for the radial component of equation of 
motion (EOM) follows from equation (93) of Paper I, but the 
second term on the right-hand side (RHS) of the latter has been 
shifted to the left-hand side (LHS). The derivation may become 
evident, when one is referred to the identity 




D 



-T-lnL, 
dr 



(20) 



In equation (11), i>k denotes the Kepler velocity (GM/r) 1 / 2 
with G and M being the gravitational constant and the mass 
of a central object, respectively. The expression (13) for the 
azimuthal component of EOM has been quoted from equation 
(102) of Paper I. 

We often refer to the expansion parameter eo = 
{v r /vip)ypF = tnat becomes very small at large radii, 
where v r and v v are the radial and azimuthal components of 
the velocity, respectively. The suffix VPF means that the quan- 
tity is evaluated by the lowest-order solution, the vanishing 
Poynting-flux (VPF) solution (see Paper I for details). 

3. Asymptotic Solutions at Large Radii 

3.1. Removal of Extended-RIAF Condition 

The original RIAF condition, 

1 dp 
a = -- — 
p dr 



GM r dp I 2 

—5- = -- — / v K = const., 

r^ p dr I 



(21) 



has been assumed in the asymptotic region at large radii. This 
condition means that the ratio of the pressure gradient force to 
the gravity should be a constant in this region, which guaran- 
tees a characteristic nature of optically thin RIAFs, i.e., a virial- 
like temperature. Combined with the equation (68) of Paper I, 
which is the lowest-order version of the r-component of EOM 



in the power series of eo, we have derived a lowest-order so- 
lution (Kaburaki 2000) in the asymptotic region. Since the 
Poynting flux vanishes identically in this solution (i.e., the VPF 
solution), it cannot explain the jet launching which is com- 
monly expected for the accretion disks of the resistive-RIAF 
type. This is because the electrodynamic launching surely re- 
quires the supply of jet-driving power through the Poynting 
flux. 

In order to overcome this difficulty, we tried in Paper I to 
improve the accuracy of the solution by taking the first-order 
corrections in eo into account. In this connection, we have re- 
placed the original RIAF condition by the extended RIAF con- 
dition, 




1(1 -n). 



(22) 



Here, n is a constant which controls the radial profiles of rele- 
vant physical quantities in the VPF solution, and in this sense 
plays a similar role to the polytropic index that replaces the en- 
ergy transfer equation (see Paper I). In contrast to the original 
RIAF condition, the extended one does not have firm grounds 
to stand on, except that it becomes identical with the original 
one in the VPF limit. Therefore, we cannot reject the criti- 
cism that it is only a makeshift policy, and if possible, such an 
obscure postulate should be avoided in deriving higher-order 
solutions. 

When the above postulate has been removed, the only re- 
maining requirement is that any improved solution should coin- 
cide with the VPF solution in the limit of vanishing corrections. 
Under such circumstances, the new policy we adopt here is to 
portion out 5ft t and D from equation (13) so as to reproduce a 
term that is proportional to v\ or v 2 , as the leading contribu- 
tion from the partial-pressure gradient force appearing in the 
identity (20). This requirement is equivalent to the condition 
(21) or (22), respectively, as far as the leading-order terms are 
concerned. Even if we follow this policy, the solution is not de- 
termined uniquely, and anyway the process of finding solutions 
becomes necessarily a kind of trial and error. 

In the following subsections, we show four examples of suc- 
cessful trials. They are all different from the outer asymptotic 
solution obtained in Paper I, which has been derived under the 
restriction of extended RIAF condition. Although only one of 
them leads to a final global solution, we dare to mention all of 
the new results. We believe that such a description will be help- 
ful for the reader to become familiar with subtle insight into the 
strategies in finding solutions and to experience how the type 
of an accretion flow (i.e., sub- or trans-critical infall) is actu- 
ally determined. If that had been omitted, the description of 
this paper would become very abrupt and less understandable. 

3.2. Following Paper I 

In specifying the radial profiles of the velocity and magnetic 
fields, it is natural to first follow the results obtained in Paper I. 
The velocity components have been written as 
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V r (0=Vin r 1/2 



D -2AC 



1/2 



2A 



(24) 



where A is a positive constant whose value is specified in the 
course of discussion, and 



2n + 1 



1/2 



GM 



(25) 



The profile of the r-component for the magnetic field has been 
fixed as 



b r (0 = B in £ 



-(3/2-n)„-(2n+l)^- 



(26) 



At the disk outer edge £ out , the magnitude of this component 
becomes comparable to the externally imposed uniform field, 
B a , as guaranteed by the relation B ln = B £out ~™- 
The above specifications result in the derivatives 



d W - 
r— m[rv 



— ln(r 2 6 r ) = 



2AcT 



1 



1 + 2AC 



dr x " 2 
and further from equations (1), (8) and (13), respectively, 



be(0 



2n + l 



B in (l + 2AC 1 )x 

x ^-(3/2-n) e -(2n+l)A«- 



SR p =^ 2 (l + 2^r 1 ) 



1/2 



and 

3?t = 



2L> 
2n+l 



'l-e- 2 ^" 1 ' 
2A 



2.4 



-1/2 



-AV 



(27) 
(28) 

(29) 



It is evident that the above set of selections for v v and b r is 
effective in keeping the expression (29) rather simple. 

Case 1. When we portion out D and 3? t from equation (29) 

as 

2n + 1 

£>= — — = const. (1/4 <£><!), 



fit 



-1/2 



,-AC 



2A 

we obtain the relation 

sft 2 2 



(30) 
(31) 



-i e -2Ar 



D 1 2n+l 111 ' 2n- 
reproducing one of our aimed form. Since the above deter- 
mined D is a constant smaller than unity for any value of n in 
the allowed range, — 1/4 < n < 1/2 (see Kaburaki 2001), the 
accretion flow in this case is said to be a 'sub-critical' 2 infall. 
Although Eg = as confirmed from 9f? t =v v> /v r , the solution 
obtained in Case 1 does not belong to a VPF solution since 
K p ^ 1 (i.e., £ r ^ 0). 



2 We follow here the terminology of the traditional accretion and wind the- 
ories (e.g., Bondi 1952; Parker 1960) in which there is a critical point in 
the governing differential equation at D = 1 (i.e., v r = Va)- Actually in 
our treatment, however, no such criticality exists at D = 1. 



The form of bm follows from the definition of 3?+ 



&„(0 = 3it& 



which yields 

r-^-lnL : 
dr 




£-(3/2-n) x 

-1/2 



-2 At 



-2(n+l)AC 



2A 



(32) 



-(1-n) 



-2(n+l)AC L -- 



i i - (i + 2,4c;- 1 ) 



-2AC 



(33) 



At large radii, the first term on the RHS of equation (33) 
is the quantity of 0(1) while the second and third terms 
are of 0(£ -1 ). They are referred to, respectively, as 

{r(d\nb v /dr)} M and {r(d\nb v /dr)} R . 

Then, shifting only the leading-order term resulting from 
(33) to the RHS in the r-component of EOM (11), we have 



SP2 

KHS = v'-vi-irj-\nb v } -± i 2 r 



dr 



M 



D 



A\ 



tinC 1 (l 



-2AC 



6A 



2n + l 

On the other hand, the LHS of EOM (11) becomes 
LHS = 



(34) 




-2 



(35) 



In the final line of the this equation, only the lowest-order terms 
in £ _1 have been kept in the curly brackets in order to match 
the form of the RHS. Finally, equating both sides, we can fix 
the value of A as 

2n + 9 

A=— — . (36) 

8(n + 2) 

Case 2. The other option for portioning out of D and 3? t is 

2n+ 1 



D 



= 2Af" 



-2AC 



-1/2 



2A 



(37) 
(38) 



in which we have b,„ of the form 



-(3/2-n) , 



-2A£~ 



-1/2 



-2nAC 
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(39) 



and the relations 

— vz = 

D 2n + 1 



K ' 



No. 



r— ln& 

dr 
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_! 1 l-(l + 2Ar 1 )e" 2Arl 



In order to keep equation (13) simple, we are naturally led 
to select 



1_ e-2.tr 1 



Differently from Casel, the accretion flow in this case is a 
'trans-critical' (see the footnote 2) infall because it starts at a 
subcritical velocity (i.e., D < 1) at the disk outer edge £ out 
(^> 1) and reaches a super-critical value (i.e., D > 1) at small 
radii (£ < 1). 

After shifting only the term including {r(dln6 v /(ir)}M to 
the RHS of equation (1 1), we have 



RHS = v I - v 
2n + l 



2 , 2(1 -n) 2 ! 
K+ 2n + l ^ 



v£ in r 1 (i 



-2A & 



-2 



The remaining LHS becomes 

2n + 9 _ . 2n-l 



LHS 



2A 



(40) 



(41) 



2(2n+l) ~ 2n + l 

according to the same procedure as in Case 1 . Equating both 
sides of this equation, we obtain 

2n + 9 



A = 



(42) 



16n ' 

which requires n > since A should be positive definite. 

3.3. Other Possibilities 

In this subsection, we discuss fairly different types of ex- 
pressions for the radial profiles of the velocity and magnetic 
fields. The new guideline in specifying them is to pay special 
attention to the identity 



4 lnF=s > or % =rlFS > 



where 
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(44) 



(45) 



These are monotonic functions as shown in figure 1 , and their 
behaviour at large radii is 

J?(0->C(l-2^- 1 ), S(0->1 (asC^oo) (46) 

within the accuracy to the first order in 
First, the radial profile of b r is specified as 

b r (0 = B in C 2 F (2n+1)/2 . (47) 

It is easy to confirm that equations (47) and (26) are equivalent 
within the accuracy to the first order in Then, we have 

2n+ 1 



B in r 2 F (2n+l)/2 5j 



(48) 



from the relation 

..2 



2n- 



(49) 



which gives 

For the radial velocity component, it turns out after some trials 
and errors that the specification 



l2n 4- 1 

v r (0=v i nr 1 = ^/-^-v K;in r 1 , 



and hence 



v 



F l/2 



(50) 



(51) 



is very interesting to examine. Then, the velocity ratio remains 
the same as in subsection 3.2, and we find very simple expres- 
sions for the magnetic Reynolds numbers: 

2D 



9? D = S, and 3? t 



F 1/2 . 



(52) 



2n+l 

Similarly to the discussion in the previous subsection, there 
are again two possibilities in portioning out D and 3? t from the 
latter of equation (52). 

Case 3. The first possibility is the specifications, 



D = 



2n+l 



= const. (1/4<D<1), 



(53) 



(54) 



which describe a subcritical accretion flow analogously to Case 
1 . In this case, we have the results 



2 ipn+1 



(55) 



and 



3ft t 2 _ 2 
— v 
D r 



2n + l 



<c 2 f 



2n- 



r—\nK = -2 +(n+l)S. 
dr 

In the last equation derived above, the two terms on the RHS of 
the logarithmic derivative of b v are both of order unity at large 
radii, since S — > 1 as £ — > oo. 

Therefore, we shift the whole term containing r(d\nb ip /dr) 
on the LHS of EOM (1 1) to the RHS: 



RHS 



~9 



d 1 S \ ^ ~2 



6A 



■{2n + 5-2(n + l)S'}r 1 ^ 



(56) 



2n + l 

where the last expression is the limiting form at large radii cal- 
culated with the aid of equation (46). Note that all the leading- 
order terms in the square brackets cancel out in this limit and 



6 



O. Kaburaki 



[Vol. , 



we have only a first order term (oc £ 1 ). The terms remaining 
on the left is 



LHS 



d 1 d I b r 
r— mv r H — r— In — 



dr 



2D dr 



2(2n + 
2n + 9 
2(2n + l) 



--{4n+10-(2n+l)S} v 2 r 



(57) 



where the last expression is also the approximate form at large 
radii. Therefore, this equation holds asymptotically at large 
radii, as far as 

Case 4. The last option in our consideration is the profiles 




(59) 



(60) 



which describe a trans-critical accretion flow analogously to 
Case 2. In this case, we have the results 



(61) 



and 



— v = 

D r 2n + l 



2 t~\ * 2 



r-^-lnL 

dr 



-2 + nS. 



As in Case 3, both terms resulting from the above logarith- 
mic derivative of b v are of order unity, and hence the term 
containing this factor in the r-component of EOM should be 
shifted altogether to the right. Then, the RHS becomes 

RHS = -1 1 { 1 + 2nS - (2n + l)^ l F) (62) 

~ -2i4 v 2 r , 

where the last expression is the limiting form at large radii. 
Again, note that the leading-order terms in the curly brackets 
have been completely cancelled out in this limit. On the other 
hand, the LHS becomes 



LHS 



VhjT* [2(2n+l) 



l -(2n + l)S}C 1 F] 



(63) 



2n + 9 



vl 



2(2n + l) r 

Therefore, the equation holds in the outer asymptotic region, 
when 

2?i + 9 



.4 



4(2ra+l)' 



(64) 



4. Global Solution 

In the previous section, we have obtained four different sets 
of asymptotic solutions at large radii. Although these sets have 
different expressions for any one of the relevant physical quan- 
tities (and their components), they are equivalent within the 
accuracy to the first order in eo (= ^ _1 ). As far as we remain 
only in this outer asymptotic region, we cannot therefore judge 
which type of the accretion flows (e.g., the sub-critical or trans- 
critical type), is more likely to fit for the reality. For the resolu- 
tion of this problem, considerations from a global point of view 
are needed. Thus, we are led to examine the behavior of the 
above sets in the opposite limit of small radius. Fortunately, as 
shown below, there is one and only one case (Case 4) in which 
the same set also serves as the asymptotic solution at small 
radii. This means that this set can be regarded as a global solu- 
tion, though the accuracy may be somewhat poor in the middle 
region. 

In order to discuss the small radius limit, we note here that 



-4-0, F^O, S^2AC\ (as £->•()). (65) 



In the first three cases discussed in the pervious section, 
the asymptotic behavior is different on both sides of the r- 
component of EOM. Indeed we obtain, LHSoc £ -3 and RHSoc 
i~ 2 in Case 1; LHSoc £~ 3 and RHSoc £~ 2 in Case 2; LHSoc 
£ -3 and RHSoc £ _1 in Case 3. However, in Case 4, equations 
(63) and (62) yield, respectively, 

271 + 1 n n 

—5— Vh n C 2 , (66) 



LHS 



RHS 



K. 



Then equating both sides, we obtain 
2n + 1 



A = 



■in 



(67) 



(68) 



which specifies the value of A in the asymptotic solution at 
small radii. 

If the values of A in the two asymptotic regions at large and 
small radii (i.e., equations [64] and [68]) coincide, the asymp- 
totic solutions match smoothly and become a global solution. 
This actually happens when 
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(69) 



2(5-Vl7) 2 

The presence of a select value of n may suggest that a prefer- 
able thermodynamic circumstance is required for realization of 
the state described by our global solution. 

The global solution indicates that the infall has a trans- 
critical nature as seen from equation (59). However, this should 
not be interpreted as a restriction on the radial profiles of in- 
fall velocity, since it is fixed always by equation (50) for both 
Cases 3 and 4. Rather, it should be interpreted as a restriction 
on the profile of the characteristic velocity Va = b r / (Anp) 1 / 2 , 
and hence, on those of b r and p. 



Other physical quantities (than v r , 



b r , be, and b v ) are 



derived straightforwardly as follows. It should be noted that all 
quantities are written in closed forms, i.e., not in the forms of 
truncated power series in We obtain from equation (14) 
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Fig. 1. Behavior of the radial parts (v r , Vg, Vtp) of the fluid velocity 
are shown together with that of the functions F(£) and S(£). The 
velocities are normalized to v± n . Although vg exceeds v r in the region 
£ $ 2, remember that the total expression for vg(£, rf) contains a small 
factor A. 



Db r 



3B- r 



2n 



and further substituting it in equation (10), 



2n + l 



(70) 



(71) 



Equation (71) indicates that the generation of a wind from an 
accretion disk is directly controlled by the parameter n (dif- 
ferently from the result in Paper 1). The direction of the wind 
is upward (i.e., vertically outgoing) when vg > 0, and down- 
ward (i.e., vertically ingoing) when vg < 0. The above selected 
value, n ~ 0.25, means the presence of a medium-strength up- 
ward wind. It can be seen in figure 1 that, as far as the profiles 
of the velocity components are concerned, they are essentially 
the same as in Paper I. 

The pressure and temperature are calculated respectively 
from equations (12) and (15) as 



p(0 



K 

8tt 



i + K t - 



8tt 



r 2 ^ 2n (i+r 2 ^) 



P 1 V K,in & -\ 



w Kp 3 K 



i+c 2 f) 



(72) 



(73) 



The deviation of temperature from the virial form is expected 
only in a middle region and remains to be rather small. 

The results for every component of the current density and 
the electric field follow from Ampere's law and Ohm's law: 

>(0-^r 2 F", (74) 



:m - 
MO -- 

E r (0 



cnB m 
4.7rr in 

C-Bin 



£-3 F (2«+l)/2 ; 



(75) 



(76) 



2n + l 2n+l V K ,mB h 



Fig. 2. Behavior of the sign-reversed radial functions of the Poynting 
flux (i.e., —P r , —Pg, —Pip). They are normalized to Vi n B? n /4:TT. 
Although —P r and — P v are everywhere positive, —Pg is positive only 
where £ < 2 and slightly negative where £ > 2. In order to see this fact 
more clearly, the formula enclosed by the curly brackets in equation 
(81) is shown as Y (£). 



(77) 



Ee(0 = - 



2n + l V KM B in C 2 F n^_ r i F ^ (?8) 



2n+l 2n + l V K ,i a B b 



4 V 3 c 

xr 3 ^ (2n+1)/2 (s-i). 



(79) 



Similarly to the case of vg, the coefficient of jg is proportional 
to n. It seems rather natural from the viewpoint of current clo- 
sure that the value of n determined in equation (69) is no-zero. 

Finally, we cite the component expressions for the Poynting 
flux, including their ^-dependences since they have been 
dropped in Paper I by accident. 

Pr(£)V) = Pr{Q sech 2 ?y tanh 2 ?7, 



P r (0 = ir E <>K = - 

47T 



2n+l VkmB? 



X 



3 4tt 
xf^fl-f'F), (80) 



Pg{i,rj) = APg(0 sech 2 /; tanhry, 



Po(£) = j-[E r b v ,-E v b r 



2n+l 2n + l V K ,i n Bf r 



x 



4 V 3 4tt 
x ^F^ll-^FS-f^FiS-l)}, (81) 



PA^'V) = P<f>(0 sec h 4 ?? tanh 2 7?, 

^(0 



^ = - A /^±l^^x 
4tt V 3 4tt 



x 



r 4 -F ,(4n+1)/2 (i-r 1 -F), (82) 



As it turns out from figure 2, P r and P v are negative every- 
where, but Pg changes sign from slightly positive (outgoing) to 
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negative (incoming) as the radius decreases across £ ~ 2. In or- 
der to see this fact more clearly, the formula in the curly brack- 
ets in equation (81) is shown as Y(£) in the figure. Taking the 
r]-dependece of P r (£, rf) also into account, we can say that the 
Poynting flux in the poloidal plane flows from the wide outer 
region (£ > 2) into the narrow inner region (£ < 2), almost along 
the surface of an accretion disk. 

It should be noted that Pg here has the opposite sign to that 
obtained in Paper I, in the inner region where Pg becomes very 
large. This means that the extrapolation of an outer solution 
that is not the global solution can actually lead to an erroneous 
result. 

5. Local Energy Budgets 

5.1. Mass Accretion Rate 

Before proceeding to energy budgets, it is convenient to in- 
troduce the mass accretion rate through a vertical cross-section 
of arbitrary radius, 



M(0 



2jrpv r r Adr] ~ 2%Apv r r' 1 



sech 4 ?7 dr\ 



2n- 



1 AB&rg, p2n 



V K , 



Mo G^F 



2/i 



(83) 



where the approximation v r = — « r sech r\ has been used, ex- 
trapolating the functional form near the equator even to large r\ 
regions (i.e., neglecting the tanh 2 r\ term). This is because the 
^-dependences are reliable only near the disk midplane owing 
to the adopted method of approximation (see Paper I). Since 
we need a co-latitudinal integration in the above calculation, 
this may cause some worry about the accuracy of the result. 
However, the most important thing is a finiteness of the inte- 
gral and its precise value does not matter on the essence of the 
following discussions. Reflecting the presence of the vertical 
flows, M varies with £ like F 2n . It is a constant only when 
n = 0, i.e., there is no wind. 

The accretion rate at the outer edge of an accretion disk is 
given explicitly as 



M = M{Ut) 



2n- 



VK.in 



(84) 



because, at this radius, F can be approximated by its asymp- 
totic form, F ~ £. First substituting B ln = Bq £^ 2 n , Vk.in = 
(GM/nn) 1 / 2 and £ out = r ou t/r; n in the above definition, and 
then solving for Bq, we obtain 

1 VgmMq 



2n + 1 a r 



5/2 



t-1/2 
>out 



(85) 



This equation will be used to eliminate Bq in the following 
subsections. 

5.2. Thermal Energy Budget 

The heating rate due to the Joule dissipation of the electric 
current per unit volume is calculated as 

•2 

.4. 



V) 



<zj (£) sech v, 



i 



Or 3.U3. I 



aA 2 
2n+l GMMq 
16ttA rL t 



Sout S 



4 F 2n (l+£-' 2 F). (86) 



In obtaining the approximate expression in the first line of 
the above equations, we have neglected jg since j r ~ j v ~ 
0(A~ 1 ), jg ~ O(l), and also neglected a term containing 
tanh 2 rj dependence. In obtaining the expression on the third 
line, a A 2 and B 2 have been eliminated with the aids of equa- 
tion (73) in Paper I and equation (85) above, respectively. 

On the other hand, the advection cooling rate per unit vol- 
ume is written as 



<?adv(0 1 



1 



V • (hpv) 
5 p d 
2 r 2 dr 
GMMq m _ 2 „ 



r 2 v r ) 



~ ?adv 

dp 
v r — 
dr 



(£) sech 4 r^, 
5 pvg 
2 r 



16ttA 



^OU 



xf¥" 



{ 1 + AnS + (An - 3)€~ 2 FS} , (87) 

where h=(5/2)(p/p) = (5/2)(p/p)is the enthalpy for an ideal 
gas. Similarly to the case of qf , the approximate expression on 
the first line of the above equations is obtained by neglecting 
the tanh 2 r\ terms. 

However, the heating and cooling terms given above does 
not balance generally. This is because i) there may exist the ex- 
change of thermal energy between neighboring fluid elements 
through conduction and/or convection (these are expressed as 
non-adiabaticities), even if the radiation loss may be negligi- 
bly small, or ii) the balance cannot be achieved in principle, 
implying the break-down of the model. The amount of this 
discrepancy is 

£dv(0-#(0 



1 GMMq 



16ttA 

x [2n i 



(i— 2nt— 4 p2n 
Sout S r x 



l + £- 2 F)(2S -1)-C 2 F{3S+1)} , (88) 

when it is expressed in terms of an additional heating (or a 
cooling if it is negative). 

There are two distinct terms in the above result; the term 
which contains n and which does not. The former dominates 
over the latter in the outer asymptotic region (i.e., in the VPF 
limit). In this asymptotic region, this term has been interpreted 
as due to the non-adiabaticities (Kaburaki 2001). Our global 
solution indicates that both terms enhance in a narrow region at 
around £ ~ 0.4. Even in such inner regions, the former drives an 
upward wind (i.e., vg > 0) when there is additive heating (i.e., 
n > 0), and vice versa, as seen in equation (71). On the other 
hand, the latter term always contributes to a cooling. Although 
it is not so certain, this fact may suggests that the latter corre- 
sponds to a radiative cooling that may enhance at such an inner 
region. 

As shown in figure 3, the advection cooling becomes largely 
exceeds the Joule heating where £ < 1. The cause of this en- 
hancement is in the monotonic increase of v r and vg toward the 
center. In principle, it is impossible to balance the advective 
cooling even by supplying heat through conduction or convec- 
tion, as far as the cooling there exceeds the peak value of the 
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Joule heating. In this sense, the present model does not seem 
to hold in the region where £ < 0.8. Thus, we are consistently 
led to the conclusion that there is an inner edge of an accretion 
disk at £ - 1. 

5.3. Electromagnetic Energy Budget 

The energy equation for the electromagnetic field is known 
as the Poynting theorem: 



du i 1 

^+V-P = -E-j = - J - (jxb)-v, 

at a c 



(89) 



where u = (1/8tt)(6 2 



(l/87r)b 2 (see Paper I) is the 



electromagnetic energy density, and P is the Poynting flux. 
The two terms on the right-most side of equation (89) represent 
losses of the electromagnetic energy through the Joule dissipa- 
tion and through work done by the magnetic force on the fluid, 
respectively. 

When we concentrate our attention mainly to the region near 
the disk midplane, we can neglect tanh 2 77 terms as before. 
The contribution from du/dt vanishes in this process, and the 
Poynting equation finally reduces to the form 

W(Z,T}) = — (3 x 6) ■ v ~ W[g) sechV 



w(0 



2n+l GMMq 
16ttA f* 



c 4 ut 2n r 5 (i 



F 2n+1 S. (90) 



As the above definition of W(£, rj) includes a minus sign, it 
means the work done by fluid on the electromagnetic force 
(hence, on the field). The explicit expression of W{£) can be 
reached either by calculating the right-hand side of the expres- 
sion on the middle line or by calculating its left-hand side di- 
rectly according to its definition, confirming that the Poynting 
equation actually holds in our case. 

Since W > at any finite radius, the fluid motion is doing 
work on the electromagnetic field everywhere in the disk. As 
seen from figure 2, the vertical Poynting flux, Pg, is positive 
in the region £ > 2 while negative in £ < 2. This implies that 
W > qf where £ > 2 and W <qf where £ < 2. Namely, in the 
wide outer region (£ > 2) the fluid in the disk is doing work that 
exceeds the local dissipation q\ and the excess is gathered into 
the small inner region (£ < 2) through the Poynting flux (see 
also the discussion at the end of the previous section). Thus, 
in the inner region, energy input through the Poynting flux can 
largely exceeds local supply through the work W, and both are 
thermalized as the Joule dissipation (i.e., W — Pe/r = 

In other words, the outer main disk is driving the innermost 
region electrodynamically, suggesting that, if a more careful 
treatment of the vertical structure and flow in this region is 
introduced in the model, the launching of an MHD jet could 
be obtained. Within the present status of our accretion disk 
model, however, this ability is wasted only on the Joule dissi- 
pation there. In a sense, the new global picture of our model is 
reasonable because it closes within itself. On the other hand, 
the previous picture stated in Paper I is less persuasive, because 
it does not close unless assuming a presence of some external 




0.5 1 1.5 2 
Normalized Radius (£,) 



Fig. 3. Behaviour of qj, g adv , — Pe /r, W and Be. Note that the for- 
mer four quantities and the latter one have different dimensions, and 
normalized to (GM M /167zAr* ut ) C^ 2 ™' and V£ .J6, respec- 
tively. 



component (e.g., an MHD jet) that is driven by the large pos- 
itive (i.e., outgoing) Poynting flux mainly emanating from the 
innermost disk. 

5.4. Binding Energy 

The binding energy per unit volume of a fluid element is 
described by the Bernoulli sum (Narayan & Yi 1994, 1995), 



, , * 1 o , GM 
2 r 



5 p GM 
2~p r~ 



(91) 



where the approximations of geometrically-thin disk and that 
of respecting midplane have been used again in obtaining the 
expression on the second line. Calculating this quantity in 
terms of our global solution, we obtain 

x {{2n + l)C\l + F) + 5C 2 F-l}. (92) 
The asymptotic values for large and small £ are 

tVsbnC 1 , as £^00, 



2n + l 
6 



(93) 



vl in C 



as £ -> 0. 



The upper line of the above equation reproduces the VPF re- 
sult (Kaburaki 2001). As seen in figure 3, the global so- 
lution predicts that B C (C) remains positive everywhere (cf. 
Narayan & Yi 1994, 1995; however see also, Nakamura 1998; 
Abramowicz et al. 2000; Blandford & Begelman 1999; 
Beckert 2000; Turolla & Dullemond 2000) and becomes di- 
vergently large in the limit of £ — >• 0, even when n = 0. This fact 
suggests that the accretion flows described by our global solu- 
tion should be ejected, at least in its fraction, before it reaches 
the center. 
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6. Summary 

Summarizing the discussions in the previous section, we 
have reached the conclusions that i) the accretion state char- 
acterized under the name of resistive-RJAF does not seem to 
extend into the region £ < 1, ii) the main disk (£ > 2) is driving 
the innermost region (£ < 2) electrodynamically, and iii) in- 
falling matter always stay unbound and cannot reach the grav- 
itational center, as a whole. Judging from these evidences, the 
most probable scenario is an ejection of the infalling matter, at 
least in its fraction, at around £ ~ 1 which may be regarded as 
the inner edge of an accretion disk. Thus, it has been shown 
clearly that one of the most preferable circumstances neces- 
sary for the MHD jet launching is actually prepared within our 
resistive-RIAF model. 
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